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Abstract 

We investigate the twist-3 effects in the polarized virtual photon 
structure. The structure functions g{ and g] of polarized photon could 
be experimentally studied in the future polarized ep or e + e~ colliders. 
The leading contributions to g[ are the twist-2 effects, while another 
structure function gl, which only exists for the virtual photon target, re- 
ceives not only the twist-2 but also twist-3 contributions. We first show 
that the twist-3 effects actually exist in the box-diagram contributions 
and we extract the twist-3 part, which can also be reproduced by the 
pure QED operator product expansion. We then calculate the non-trivial 
lowest moment (n = 3) of the twist-3 contribution to gj in QCD. For 
large N c (the number of colors), the QCD analysis of the twist-3 effects 
in the flavor nonsinglet part of gj becomes tractable and we can obtain 
its moments in a compact form for all n. 
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1 Introduction 



In recent years, there has been growing interest in the study of spin structures of 
photon. Especially, the first moment of the polarized photon structure function gj 
has attracted much attention in the literature @-|| in connection with the axial 
anomaly, which is also relevant to the analysis of the nucleon spin structure function 
giiuci rj^g p j ar j zec [ photon structure functions may be extracted from the resolved 
photon process in the polarized version of the ep collider HERA || |7|. More directly, 
they can be measured by the polarized e + e~ collision experiments in the future linear 
colliders (Fig.l), where — Q 2 (— P 2 ) is the mass squared of the probe (target) photon. 

For the real photon (P 2 = 0) target, there exists only one spin- dependent struc- 
ture function, gJ(x,Q 2 ), which is equivalent to the structure function W±(x, Q 2 ) 
(gj = 2W2) discussed some time ago in || [J. The leading order (LO) QCD correc- 



tions to g{ for real photon target was first calculated by one of the authors in [1C 



and later in |TT, 0], while the next-to-leading order (NLO) QCD analysis has been 
performed in |T2|, [TB[ . 

In the case of virtual photon target (P 2 ^ 0), and especially when we deal with 
the kinematical region, A 2 C P 2 < Q 2 , where A is the QCD scale parameter, we 
can calculate the whole structure functions up to NLO in QCD by the perturbative 
method, in contrast to the case of real photon target where in NLO there exist non- 
perturbative pieces [|14], |15|. In this context, the spin-independent structure functions 
F^ix, Q 2 , P 2 ) and F2(x, Q 2 , P 2 ) as well as the parton contents of unpolarized virtual 
photon were studied in the above kinematical region in LO [|16|] and in NLO [[17]- 
The target mass effects on both the unpolarized and polarized virtual photon 
structure functions were discussed in LO |p3| . More recently, the spin-dependent 
structure function gj (x, Q 2 , P 2 ) of virtual photon has been investigated up to NLO 
by the present authors in [p4| , and also in the second paper of |13| . 



Generally, for the virtual photon target, there exists another structure function 
g\ (x, Q 2 , P 2 ), which is the analogue to the spin-dependent nucleon structure function 
g2 Ucl . In the language of operator product expansion (OPE), it is well known 
that both twist-2 and twist-3 operators contribute to g™ 01 m the leading order of 
l/Q 2 - The same is also true for gj. In this paper, we shall investigate the twist-3 
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contribution to gl in the leading order in QCD, and show that they are sizable in 
contrast to the nucleon case, where the experimental data show that the twist-3 
contribution appears to be small [35, [26| . 

In the next section we discuss the two structure functions, gj and gl, which 
describe the deep inelastic scattering off polarized virtual photon target. (How the 
information on gl and gl can be extracted from the spin dependent experiments of 
e + e~ — > e + e" + hadrons is shown in Appendix A.) They are related to the s-channel 
helicity amplitudes which appear in the forward virtual photon-photon scattering. 
We examine the box diagrams for the photon-photon scattering to obtain the parton 
model predictions for gl and gl- We then extract a piece which is a deviation from 



the Wandzura-Wilczek relation |27j. In Sec. 3 we show in the framework of the 
pure QED operator product expansion that the extracted piece actually arises from 
the twist-3 effects. In Sec. 4, we examine the QCD twist-3 effects in gl in LO. We 
calculate the non-trivial lowest moment (n = 3) of the twist-3 contribution to gl. 
For large N c , the QCD analysis of the twist-3 effects in the flavor nonsinglet part of 
gl becomes simple and we can obtain its moments in a compact form for all n. The 
final section is devoted to the conclusion and discussion. 



2 gl (x, Q 2 , P 2 ) and box-diagram calculation 

Let us consider the virtual photon-photon forward scattering for 7(g) + j(p) — > 
7(9) + lip) illustrated in Fig. 2, 

T^ P r(p,q) = i J d 4 xd%d^ x e^y- z \0\T(J^x)JMUy)Mz)M- (2.1) 

where J is the electromagnetic current, and q and p are four- momenta of the probe 
and target photon, respectively. Its absorptive part is related to the structure tensor 
W^ pT (p, q) for the photon with mass squared p 2 = —P 2 probed by the photon with 
q 2 = -Q 2 : i 

WWO,*?) = -lmT„ upT (p,q) , (2.2) 
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The anti-symmetric part, W^ upr , under the interchange // <-> v and p «-> r, can be 
written in terms of the two spin-dependent structure functions, gj and g\ as 



Kv P r = t^K'-W gl - (J-W s 2 7 ], (2-3) 

where the two tensors I_ and J_ are explicitly given by 

(I-)pupr = V Q^nuXa^pT^q^p, (2-4) 

{J-)i±v P t = e^x^prapq x p a q a p (} -p- qt l Mu\a£pT a(5 q x pp- (2.5) 

In fact we observe that 7_ and J_ are related to those of eight independent kinematic- 
singularity-free tensors, introduced by Brown and Muzinich (see Eqs.(A3)-(AlO) 
of Ref.|§) to express the virtual photon-photon forward scattering amplitude, as 
follows: /_ = I 2 — h and J_ = I 7 — I s . 

It may be useful here to see the relations between structure functions gj (i= 1,2) 
and s-channel helicity amplitudes, which are defined as 

W{ab\a'b') = el{a)e* p {b)W^e v {a')e T {b') , (2.6) 

where e^(a) represents the photon polarization vector with helicity a, and a, a' = 
0, ±1, and b, b' = 0, ±1. They are related as follows: 

/ 2 2\l/2 

{iy(ii|ii) - w(i — 1 1 1 -1)} - p g {^(ii|oo) + w{oi\ -10)} 

{W(11\U) - W(l -1|1 -1)} - ^^{W(11|00) + W(0l\ -10)} 

(2.7) 

where X = (p ■ q) 2 — p 2 q 2 . 

The photon structure functions gj and gl are just the analogues to the nucleon 
counterparts g± ncl and g 2 ncl , respectively. But it is noted that g 2 exists only for the 
off-shell or virtual photon (P 2 7^ 0) target. 



9i = 


1 

2X 


gl = 


-1 


2X 
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Now let us calculate gj and gl in the simple parton model by evaluating the box 

fiupr 



diagrams depicted in Fig. 3. We introduce two projectors (P/) / " /pT and [Pj)^ vpT 



{PiY vpT 
{PjY vpt 



AX 2 
1 

4X2 



1 + 



p 2 q 2 
2{p-qf 



y v p t + -(j Y upT 



jlVpT 



+ 1 + 



(p-q) 2 

2p 2 q 2 



(J- 



(2.8) 
(2.9) 



which satisfy 



P/.(/_) = l, Pj.(J_ 
Pj ■ (/_) = 0, P, • (J_ 



1 . 



(2.10) 



Then g± and #2 are given by 
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{p-q)\PiT pT W A IMvpT , 

-ip ■ q) 2 (p J r pT w A ^ 



(2.11) 



Applying the projectors Pi and Pj to the box-diagram (massless quark-loop) con- 
tributions, and ignoring the power corrections of P 2 /Q 2 , we obtain 



7(Box) 



9i 



^(Box), 



x,Q 2 ,P 2 



3a 



71 



N f (e 4 ) 



q 2 
p 2 



(2x - 1) In ^ - 2(2a; - l)(lnx + 1) 



(2.12) 



3a 



x iQ ->P ) — Nf(e ) 

71 



o 2 

(2x - 1) In ^- + 2{2x - 1) In a: + 6x - 4 
pi 



(2.13) 



where x = Q 2 /{2p ■ q), a = e 2 /47r, the QED coupling constant, and (e 4 ) 



Y^fJi^f/Nf, with Nf being the number of active quark flavors. The above results 
are consistent with those of Ref. [28], where polarized gluon structure functions were 
considered. It is noted that the first moment of #2 vanishes, i.e., #2 satisfies 



the Burkhardt-Cottingham (BC) sum rule [29]: 



/ dxg 
Jo 



7 (Box) 



x,Q 2 ,P 2 ) = 0. 



(2.14) 



We will see from the OPE analysis in the next section that the BC sum rule for 
g\ generally holds in the deep-inelastic region, Q 2 ^> P 2 . Also note that the sum 



9i 



7 (Box) 7 (Box) 



+ 92 does not have lnQ 2 /P 2 behavior. 
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In the case of nucleoli, the spin-dependent structure function g^ 1 receives both 
twist-2 and twist-3 contributions, 



(x, Q 2 ) + 



(x,Q 2 ) 



(2.15) 



and the twist-2 part of g% is expressed in terms of g™ by so-called Wandzura- 
Wilczek (WW) relation |7|: 

micl.tw.2 / 



92 



(x,Q 2 



gT cl ' WW (x,Q' 



= -#i ucl 0, Q 2 ) + f 1 —9i ncl (y, Q 

J x y 



(2.16) 



Thus the difference, g™ cl = g™ 01 — <?2 Ud ' > contains the twist-3 contribution only. 
The experimental data so far obtained show that the twist-3 contribution to g™ ' 



appear to be negligibly small |25], |26[ . 

Now we may ask what about the photon structure function g\: Does g{ also 
receive twist-3 contribution? If so, are they small like the nucleon case, or, sizable? 
Does the WW relation also hold for gj, in other words, is the twist-2 part of g% 
expressible in terms of gjl These issues will be discussed in the next section. 

Here let us apply the WW relation to the results of box-diagram calculation, 



gJ (Box) and g 2 7(Box) in Eqs.(^T2}gT3|), and define 



7WW,(Box) 
92 1 



x,Q 2 ,P 2 



_^7(Box) ( 



x,Q 2 ,P 2 ) + 



dy 

x y 



gl iBo *\y,Q 2 ,P 2 



[2.17) 



Then we find that the difference, g^ Bca ^ = gj ~ ^^ ww '( Box ^ i s given by 



_7(Box) = ^ jV/(e 4 ) 
7T 



o 2 

(2a; - 2 - \nx) In - 2(2x - 1) lnx + 2(x - 1) + ln 2 a; 
pi 



Its n-th moment, g^n = Jo dxx n l g1^ B {x, Q 2 , P 2 ), is 

1 . Q 2 2 



^ ^ Nf{ei) n-l 



71 



71 



n(n + l) P 2 (n + l) 



n- 



(2.18) 



(2.19) 



In Fig. 4, we have shown the Box diagram contributions to the structure functions, 
(?7(Box) , g\ ( ' Box ' ) as well as the gl^ ^ given in the above equation (|2.18| ) as functions 
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of x for Q 2 = 30 GeV 2 and P 2 = 1 GeV 2 . We can see that g^"^ is comparable 
in magnitude with g^ Bov ^ for large region of x. Now it is well expected by analogy 
with the nucleon case that g^ 3 ^ arises from the twist-3 effects. In the next section 
we will be convinced that g^"^ is indeed the twist-3 contribution. 

3 OPE analysis and pure QED effects 

Applying the OPE for the product of two electromagnetic currents, we get for the 
fi-u antisymmetric part 

i J d 4 xe^T(J,(x)J u (0)) A = -te^q x ]T ■ ■ ■ 

(3.1) 

where i?" 2 ^ and are the twist- 2 and twist-3 operators, respectively, and and 
Efyi are corresponding coefficient functions. The twist-2 operators R^ 2 )i have totally 
symmetric Lorentz indices cr/ii • • -// n _i, while the indices of twist-3 operators i?" 3 ^ 
are totally symmetric among fj,i • • • fj, n -i but antisymmetric under a <-> /ij. Thus 
the "matrix elements" of operators and i?" 3 ^ sandwiched by two photon states 
with momentum p have the following forms: 

(0\T(A p (-p)R^r^A T (p))\0) Amp = -m« 2)i e pra V 1 • " -jA-V 

— (traces) , (3.2) 
(0|T(^(-p) J R^-^- 1 A r (p))|0) Amp = -^ 3)l e pr > 1 • • .p"»- V 

— (traces) , (3.3) 

where the suffix 'Amp' stands for the amputation of the external photon lines and 

i _ n— 1 

{"■^i . . . nMn-i} _ 



-pro? P • • ■ p 

n 



WV 1 • • "P^- 1 + E • • -P" • • •P' 1 - 1 ] >(3-4) 

.7 = 1 



n-1 ., 1"- J 



-r [f> „{M1 ] . . . „Mn-l} -r Ml „Mn-l \" r Wl,.,/.,. „Mn-l 

(3.5) 



n n j=1 
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Using Eqs.( |3.2| )-( |3.5| ), we can write down the moment sum rules for gj and g\ as 



dxx^gJ^Q 2 , P 2 ) = £af 2> i^ 2)i (Q 2 ) 

i 

n — 1 



dxx^gKx, Q\ P 2 



n 



E«(2)^(2)«W 2 )+E«(3)^r3),(Q 2 ) 



(3.6) 

(3.7) 



From this general OPE analysis we conclude: 

(i) The BC sum rule [T2§] holds for gj, 

[ 1 dxg2(x,Q 2 ,P 2 ) = 0, 
Jo 

(ii) The twist-2 contribution to g% is expressed by the WW relation 
" 1 E^hi^m = fdxx^g^^Q 2 ^ 2 ) 



n 



with 



7 ww 

92 



>, Q 2 , P 2 ) = -gj(x, Q 2 , P 2 ) + ^ %J(y, Q 2 , P 

Jx y 



(iii) The difference, g\ = g{ — g\ WW , contains only the twist-3 contribution 

n — 1 



Jo 



dxx n - x gl{x,Q\P 2 



n 



(3.8) 

(3.9) 
(3.10) 

(3.11) 



Let us now analyze the twist-3 part of g\ in pure QED, i.e., switching off the 
quark-gluon coupling, in the framework of OPE and the renormalization group (RG) 
method. In this case the relevant twist-3 operators are the quark and photon oper- 
ators, which are given, respectively, by 



j-,(T/il---/i n _l 

"(3)7 



i n ~ x e\ ^757 [ct '/J {mi] • ■ ■ D^-^ip - traces , (3.12) 
\i n ~ x e [(T ' aM F a{ ^d^ ■ ..d^' l} F^< - traces , (3.13) 



where e q is the quark charge, = d^+ieA^ is the covariant derivative and { } means 
complete symmetrization over the indices, while [a, fij] denotes anti-symmetrization 
on afij. With the above photon operator R%\ 7 , we have a? 3 ^ = 1 . The coefficient 
functions corresponding to operators R%\ q and -R™ 3 ) 7 , 



(3.14) 



7 



satisfy the following RG equation to lowest order in a, 







-pn ( Q 2 



a 



7, QED (« 



Q 1 



(3.15) 



where 7^ ED (a) is the anomalous-dimension matrix. To lowest order in a, this matrix 
has the formP] 

(3.16) 



7„ QED («) 





'477-^(3)9 



Here K? 3 \ q represents the mixing between the photon operator anc ^ the Q uar k 
operator R%\ q - Evaluating the triangular diagrams given in Fig. 5 and taking into 
account the color degrees of freedom in the quark-loop, we find 

1 



The solution of ( |3.15|) is 

^(3)(^,« 



A (3k 



cxp 



-2Ae A — 

q n(n + 1) 



-7 QED («)ln^ 



^3)(l, 



a 



(3.17) 



(3.18) 



To lowest order in a, the exponential and the coefficient functions E ^(l,aj are 
written, respectively, as 



exp 



-i 7 QED (a)ln^ 



(3.19) 



a) 



l + 0(a) 



E, 



(3)7 



1, a) 



£ 3e ^"3)7 • 



(3.20) 



The "matrix element" a%\ q of the quark operator R%\ q between the photon states 
is obtained by evaluating again the triangular diagrams in Fig. 5 and expressed as 



a 



1 



P 2 



^[-l K ^Jn— + 3etA n 



Air 



'(3)g 



q^-{3)q 



Inserting Eqs.( p.l7| )-( |3.21| ) into ( |3.11| ) and remembering 
n-th moment of ~g\ lIy pure QED, 



(3.21) 

1, we obtain for the 



7, n I 

92 I QED 



n - 1 a 4 
n An q 



A Q 2 

■^TT) ln p- 2+A ^ + B ^\ 



(3.22) 



1 We follow the convention used by Bardeen and Buras |1J] to write the mixing anomalous 
dimensions between the photon and other operators. 
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The dependence on the renormalization point fi disappears. And we note that 
although ^4(3)q and B™ 3 ^ are individually renormalization-scheme dependent, the 
sum ^4(3)^ + -S(3) 7 is n °t [0- The calculation of box diagrams in Fig. 3 gives 

A ^ + = 8 {(^TTp " h) ■ (3 - 23) 



Now adding all the quark contributions of active flavors and replacing 3e^ in ( |3.22| ) 
with 3iV/(e 4 ), we find that the result is nothing but <72n given in (|2.19|) which is 



derived from the box-diagram calculation. Thus it is now clear that <?2n 18 indeed 
the twist-3 contribution. 



4 QCD effects 

We now switch on the quark-gluon coupling and consider the QCD effects on the 
twist-3 part of g\. In the case of nucleon, the analysis of #2 j the twist-3 part of the 
structure function g^ 1 , turns out to be very complicated pl|-[]43|. This is due to 
the fact that the number of participating twist-3 operators grows with spin (moment 
of ^2 Ucl ) an d that these operators mix among themselves through renormalization. 
Therefore, the Q 2 evolution equation for the moments of cannot be written in 
a simple form, but in a sum of terms, the number of which increases with spin. The 
same is true for ~g\- 

Writing down the coefficient functions E^J^p, g 2 , an, which correspond to the 
relevant twist-3 operators R^)i contributing to g2, in a column vector (Q?, g 2 , a 
the RG equation for E t 3 \ can be written to lowest order in a as fli"4] 



( d d \ n Q 2 n Q 2 



a 



(4.1) 



where f3(g) is the QCD (3 function and 7„ is the anomalous dimension matrix. The 
solution is given by 



■pn (Q 2 



Texp 



m 2 ) 



dg 



,l n(g' ,a) 



ft 



(4.2) 
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The twist-3 photon operator is again given by Eq.( |3.13| ). In the convention 

we use now, where the photon coefficient function E™ 3 ^ is set at the bottom of the 
column vector E r 3 \, the matrix 7„ to lowest order in a has the form 



7n 



K n (g 2 ,a) 







(4.3) 



where j n represents the mixing among hadronic (quark and gluon) operators and a 
row vector K n describes the mixing between the photon operator R^y f and other 



hadronic operators. Then the evolution factor is given by [nj 



Texp 













(4.4) 



with 



a , K n {g t2 ,a) 



' . „ W 2 ) 

dg wo" 



(4.5) 



Expanding K n (g 2 ,a), (3{g) and ■y n (g 2 ) in powers of g, 



K n (g 2 ,a) 

Ug 2 ) 



^i?f + 0{ag 2 ) 



g' 



16tt 2 

.2 



9 —¥^+0(g% 



167T 



(4.6) 
(4.7) 
(4.8) 



we see that the dominant contributions, which behave as InQ 2 , are coming from X n , 
in other words, from the photon coefficient function E^ ry (Q 2 / /i 2 , g 2 , a) . Inserting 
the solution Efy and a™ 3 ^ 7 = 1 into Eq.( [3.11| ), we obtain the following formula for 
the n-th moment of ~g\ in LO: 



dxx n - x gl{?M\P 



n — 1 2-na 

n j3 Q 



(0)1. 



a 2 dg 



12 



(g' 2 y 



exp 



2A g a . 



E ( " 3) (l,0)]. ,(4.9) 



where i and j run over hadronic (quark and gluon) sector only. 
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The evaluation of the n-th moment of g^ is feasible when n is a small number. 
But as n gets larger, it becomes more and more difficult a task due to the increase 
of the number of participating operators and the mixing among these operators. 
However, we will see that in a certain limit the analysis of the moments becomes 
tractable. In the following subsections we consider the two cases: (1) the non-trivial 
lowest moment (n = 3) of gl] and (2) the flavor nonsinglet part of ~g\ for large N c . 
In case (1) the number of the participating operators is limited, and we can get 
all the information on the necessary anomalous dimensions. Thus we obtain the 
LO QCD prediction for the third moment of ~g\. In the QCD analysis of photon 
structure functions, the contributions are divided into two parts, the flavor singlet 
and nonsinglet parts. In case (2), we show that in the approximation of neglecting 
terms of order 0(1/ N 2 ), we can evade the problem of operator mixing for 'g~1^ NS \ 
the flavor nonsinglet part of gj, an d obtain the moments of g1^ NS ^ in a compact 
form for all n. 

4.1 The third (n = 3) moment of gl 

Let us start with the analysis of the flavor nonsinglet part. Besides the photon oper- 



ator given by Eq.( |3.13| ), the following four types of twist-3 operators contribute 
to^: 

^r-Mn-i = ^-l^ 757 ^ jD W... jD Mn-ig^_^ races ^ ( 410 ) 
Rl^'"^- 1 = ^{(Vi-Vn^i + Ui + Un^ + iVt-Vn-^i + Ui + Un^t)}, 

(Z = l,. ..,n-2) (4.11) 
R^'^- 1 = i n - 2 mS'^ 5 - f a D^---D^- 2 Y n - 1 Q ch ^- (traces), (4.12) 

,71—1 „, r ^- 

+V>(i V> - m)-f 5 j a D^ ■ ■ ■ D l *- a 'f n - 1 Q ch il>] - (traces) 



RT 1 ^- 1 = i n ~ 2 S'^^D^ ■■■D^- 2 Y n - 1 (iP -m)Q ch i) (4.13) 



with 



Q ch = Q 2 -(e 2 )l , (4.14) 

where Q is the Nf x Nf quark-charge matrix, (e 2 ) = YU^i e l/-^f anc ^ 1 ^ s an Nf x 
Nf unit matrix with Nf being the number of active flavors, and m represents the 
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quark mass. The symbol S' denotes symmetrization on the indices /X1//2 • • ■ l^n-i 



and antisymmetrization on afi%. The operators in Eg. ( OT ) contain the gluon and 
photon field strength and F^, and their dual tensors = \e^ va fiG a!3 and 
F^w = \^iLvaf3F afi . Explicitly they are given by 

y x = +i n gS^ 5 D" 1 ...G aw --- J D^-V" _1 <5 c V- (traces), (4.15) 

= -i n - l gS'^D^ ---G^ 1 ■■■D^- 2 ^ Ln - 1 Q ch i)- (traces), (4.16) 

V t = -i n eS^ 5 D^ ■ ■ ■ F an ■ ■ ■ D» n - 2 ^n-iqch^ _( traces ) ; ( 417 ) 

Ui = -\-i n ~ 1 eS'~tpD fJ ' 1 ■ ■ ■ F a ^ 1 ■ ■ ■ D^ n ~ 2r y^ n ~ 1 Q ch tp — (traces), (4.18) 

where g and e are the QCD and QED coupling constants, respectively. The operator 
K% in Eq.( ^Tj ) is proportional to the equation of motion (EOM operator) [^, M . 

We emphasize that V\ and U\ (regardless of quark charge factor Q ch ), which are 
not present in the twist-3 contribution to nucleon structure function g™ 01 ; m ust be 
included in the analysis of ~g\- The reason is that we are here considering not only 
QCD but also QED, and thus the covariant derivative should read as 

D lt = d ll -igA a ll T a + ieA lt , (4.19) 

where A® and A^ are gluon and photon fields, respectively, and T a is color matrix. 
Then the commutator, 

[Dfj,, D u ] = -igGlJ a + ieF^ , (4.20) 

leads to the appearance of VJ, Ui, Vj, and Ui terms. As far as the mixing anomalous 
dimensions among the hadronic operators, i.e., those given in (|4.10|) -( 4.13|) , are 



concerned, terms V\ and U\ are irrelevant. But they are indispensable to the correct 
evaluation of mixing anomalous dimensions K n ^\ between the hadronic operator R%v 
and photon operator R%yy We need to have K n> i of order O(o) in the leading 
logarithm approximation, but R? 3 \ t without V\ and TJ\ terms gives K„ ; i ~ (D(g 2 a), 
since VJ and Ui terms already have the QCD coupling constant g. In Appendix B we 
calculate the mixing anomalous dimensions \ of order 0(a) for arbitrary n and 
show that Ui term (but not V{) indeed plays an essential role. Another important 
consequence of introducing ieA^ into the covariant derivative is that with this 
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new term we can show that the photon matrix element of EOM operator R E , more 
precisely, (0\T(A p (— p)R^ 1 E ' fln ~ 1 A T (p))\0) Amp actually vanishes at 0(a). 

The twist-3 hadronic operators given in ( f4.10| )-( f4.13| ) satisfy the following rela- 



tion gg, |37|, |39|, |40 



n - 1 n " 2 

(3)F ~~ — " - ft (3)m Z^V i 1 t J-"'(3)Z - rt (3) J B • 

' i Z=l 

Hence, including photon operator, there are, in total, n + 1 independent operators 
which contribute to the n-th moment of gl . We have a freedom in choosing 
hadronic operators as independent bases. But we should keep in mind the following: 
due to the constraint Eq. ( |4.21| ), a different choice of operator bases assigns different 
values to the coefficient functions at the tree-level, which was first pointed out by 
Kodaira, Yasui and one of the authors |3£j. In the basis of independent operators 
which includes R%\ F but not R%\ m , the tree level coefficient functions are given by 

£ ( V(tree) = l, ^(tree) = . (4.22) 
On the other hand, if we eliminate R%\ F , we have 

77 — 1 

£ ( n 3)m (tree) = — — , Entree) = n - 1 - Z . (4.23) 

We always have Efy E (tree) = . So a different choice of the operator bases leads 
to different forms for the anomalous dimension matrix and the coefficient functions 
but the final result for the n-th moment of gl should be the same (See Appendix 
C). 

Now we take n = 3 and evaluate the third moment of ijl^^ . From now on we 
omit the superscripts n = 3. The relevant hadronic operators are four: R(3)f, R(S)i, 
R(3) m and R(3)e- Let us take Rmi, R(3) m and R(3)e as independent operators. In 
these operator bases, the tree level coefficient functions are given by 

£ (3 )m(l,0) = ! £ (3) i(l,0) = l . (4.24) 

The 3x3 anomalous dimension matrix 7^ for hadronic operators has a form, 

/Ti? \ 
7 (0) = €\ ¥± , (4.25) 
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with M 



€1 = -\C F , ^\ = - l -C F . (4.26) 



Note that we follow the convention of Bardeen and Buras |Tj] to define the anoma- 
lous dimension matrix. The matrix -y^ is triangular and, therefore, its eigenvalues 
are 7^ , 7mm 5 an d Iee- m ^ ac ^ we onr y need the information on the upper- left 2x2 
submatrix for the analysis, which is decomposed as 

7 (0) l( 2 x 2 ) =l!?Pi+€lP2, (4.27) 
where Pi and P2 are projection operators and given by 



with 

7ml —ZC t 



1?}-1 { ± 3(9CW9C F ) 



(4.29) 



The anomalous dimension K^ m is found to be null for all n. So we have K$ = 0. 

Inserting these informations into the moment formula for g% in Eq. (|4.9|) , we find 
for the third moment of the nonsinglet part g~1^ NS \ 



V$g3 = f dxx*gr s \x,Q*,P* 
Jo 



2 a 2tt 




1 l + lT S /2(3o 

<P 2 )) 



(4.30) 



where we have revived the superscript NS and k[°^ NS is obtained from Eq. (p.8|) in 
Appendix B as 

^ = _ 24iV/((e4) _ (e2)2) _l__ (431) 
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The above result for #2 n^s * s indifferent to the choice of an independent set of oper- 
ators. In Appendix C we take R(3)f,R(3)i,R{3)e as independent operators, replacing 
R(3)m with -R(3)f, and show that we obtain the same gl }. 

Let us move to the third moment of the singlet part gl . For n — 3, there are 
five independent hadronic operators contributing to ~g~2^ S \ apart from the photon 
operator R%\ 7 - Following the work of Kodaira et.al.|44[], we take R(L l7 R?3) m , R(3)ei 



T( 3 )_b, and Tr 3 ) E for an independent set. The first three are the analogues to R(g\i, 
R(3) m , R{3)E in the nonsinglet case, obtained by replacing the quark charge factor 
Q ch with an Nf x Nf unit matrix 1. The rest are the BRST-exact and the gluon 
EOM operators, respectively, whose explicit expressions are given in Ref . ||4 ] . 
The tree level coefficient functions corresponding to these operators are 



£(3)1 (1,0) 



2 
3 



(4.32) 



and others are zero. The mixing anomalous dimensions among these operators have 
been calculated and form a 5 x 5 matrix. The physically relevant part is the following 
2x2 submatrix: 

.(0)5 



with 32, 33, 44 



V 



7ff 

.(0)5 
, / ml 



A0)s 

I lm 



A(0)S 
/ mm 



(4.33) 



/ml 



6C G 
4 



9 



-C F + -N f , 
3 3 f ' 

.(0)5 = 12C 
J Imm i-^^JJr 



*(o)s 

Tim 



(4.34) 



The matrix triangular and, therefore, the same procedures as the nonsinglet 

case can be applied here. We obtain for the third moment of the singlet part g?^ , 



77 7(5) 
y2,n=3 



dxx 2 gl^ S \x, Q 7 



P 2 



2 a 



2tt 



3 4tt (3 a s (Q 2 



(0)5 



<e 2 > 



1 + 1^ S /2P 



x 



' a s {Q 2 ) \ 

MP 2 )) 



7™ S /2/3o+l' 



(4.35) 
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with 

K{ 0)S = -2AN f (e 2 )^ . (4.36) 

In Fig. 6, we have plotted the Q 2 evolution of the flavor singlet, nonsinglet 
components as well as the total of the third moment of g^{x, Q 2 , P 2 ) in units of a /it 
for P 2 = 1 GeV 2 with Nf = 3. The flavor singlet and nonsinglet components show 
somewhat different behaviors of Q 2 dependence, and the singlet component gives 
larger contribution due to the charge factor. 

4.2 Flavor nonsinglet part of for large N c 

For the case of flavor nonsinglet nucleon structure function g2 UC ^ NS \ it has been 
observed by Ali, Braun and Hiller (ABH) |45|] that in the large N c limit the twist-3 
part, g~2 Ucl ^ NS \ obeys a simple Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 



equation [46|. In their formalism of working directly with the nonlocal operator con- 



tributing to the twist-3 part of g^ nc ^ NS ^ ) they showed that local operators involving 
gluons effectively decouple from evolution equation for large N c , which is the number 
of colors. Later, the ABH result on g^ nc ^ NS ^ was reproduced by one of the authors 
[ fT?H in the framework of the standard OPE and RG method. At large N c , the op- 
erators involving gluon field strength G^ v decople from the evolution equation of 
_nuci(ivs)^ anc j wrio i e contribution in LO is represented by one type of operators. 
The same is true for the flavor nonsinglet part of g%. 

Let us take R%\ F , R&v, R%\e as independent hadronic operators, eliminating 
R%) m - The advantage of this choice of operator basis is that from Eq.( |4.22 ) the 
hadronic coefficient functions take simple forms at the tree-level [39fl, 



E ( " 3)F (1,0) = 1, E^ ){ (1,0) = forZ = l,-..,n-2. (4.37) 

The mixing anomalous dimensions for these operators are very complicated. How- 
ever it was found in Ref.|j4"7| that, in the approximation of neglecting terms of order 



0(1/ N 2 ) and thus putting 2Cp = Cq, the (F,F) and (l,F) elements are reduced 
to have simple expressions: 

11 n 1 

1 ( Sf = 8C F (S n ---—), with S n = £-' (4.38) 

4 zn . =1 j 
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7$f = 0, for Z = l,...,n-2. (4.39) 

Note that the corrections are of e>(l/iV 2 ), about 10% for QCD (N c = 3). 

Inserting the above results ( |4.37| - |4-39| ) into Eq. (fO|) , we find that, for large N c , 
the n-th moment of g1^ NS ^ is given by 



' , n-l--y(NS)^ n2 _ n ~ l a 27r ^(0) 

(0) 



'oc s {Q 



n 4tt /WQT" n, *'l+ 7 &/2A> 

2N\7i°W/2/3o+l' 



X < 



a s (P 2 



(4.40) 



with 



K® = -24iV / ((e 4 ) - (e 2 ) 2 )^^^ • (4.41) 

We now perform the Mellin transform of Eq. (|4.40|) to get jfl^ N (x, Q 2 , P 2 ) as 
a function of x. The result is plotted in Fig. 7. Comparing with the pure QED 
box-graph contribution, we find that the LO QCD effects are sizable and tend to 
suppress the structure function #2 both in the large x and small x regions, so 
that the vanishing n — I moment of yl^^, i.e. the BC sum rule, is preserved. 

As for t ne flavor singlet part of ~g~2, it is expected that a similar simplifi- 

cation may occur for large N c and its moments may be written in a compact form 
for all n as in the case of ■ At the moment we do not know how to solve the 

mixing problem in the flavor- singlet sector to get an analytically simple formula for 
the moments of <?2 for large N c . This is an interesting subject which should be 
pursued. 



5 Conclusion 

In the OPE of two electromagnetic currents, we expect the presence of the twist-3 
operators in addition to the usual twist- 2 operators. From the study of the lepton- 
nucleon polarized deep inelastic scattering, we have learned that the twist-3 contri- 
bution does not show up as a sizable effect, since the nucleon matrix elements of the 
twist-3 operators are found to be small in experiments. 
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In this paper, we have investigated the twist-3 effects in g\ for the virtual photon 
target, in the pure QED interaction as well as in the LO QCD. We have found that 
the twist-3 contribution is appreciable for the photon case in contrast to the nucleon 
case. In this sense, the virtual photon structure function g\ provides us with a good 
testing ground for studying the twist-3 effects. We expect that the future polarized 
version of the ep or e + e~ colliders may bring us important information on polarized 
photon structure. More thorough QCD analysis including the flavor-singlet part is 
now under way. 
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Appendix 



A Two-photon process e + e — > e + e + hadrons 

The information on the polarized structure functions gj and g 2 can be extracted 
from the experiments of the two-photon annhilation process with polarized e + e~ 
beams as shown in Fig. 1, 

e ± (/i)e T (/ 2 ) -> e ± (l' 1 )e T (l' 2 ) 1 (q) 1 (p) -> e ± (l' 1 )e T (l' 2 ) + hadrons, (A.l) 

with the virtual photon momenta, q — h — l[ and p = l 2 — 1' 2 . The cross section for 
this process is written as |48[| , 



P _ (4tt Q ) 2 (27r)^(p + g - P x )dT d\d% 

da ~ p 2 q2 Pi( P oi)P2( P oi)M» p M UT 4[(/i /2)2 _ m4]1/2 2E[2E 2 (27rr 1 ' 



where m is electron mass, E[ 2 are scattered electron (positron) energies, Px = ^iPi 
and dT = U i d 3 pi/2p i Q(2n) s are the total momentum and the phase-space volume, 
respectively, of the produced hadron system, and M VT is the transition amplitude for 
77 — > hadrons. For the polarized e + e~ beams, the photon density matrices PiL n 

and palpal) are s iven b y 



n^ u - 1 Tr 



= 2tme^s laqf3 /(-q 2 ) , (A.3) 
PSU = 2ime^s 2a p p /{-p 2 ) , (A.4) 

where si^ are the initial e + (e~) polarization vectors. When the incident beams are 

longitudinally polarized and at high energies, then polarization vectors are expressed 

as ^ 

*? = /k- (* = 1,2), (A.5) 
m 

with /ii = ±1, representing the helicity states of the beams. 

The absorptive part of the 77-forward scattering amplitude W^ upr in (|2.2|) is 
related to the following integrated quantity over the phase-space volume of the 
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produced hadron systemQ: 

(A7ia)W, upT (p, q) = ^J M; p M UT (2ir) 4 5(p + q - P x )dT . (A.6) 
Applying p^ pol) and p p 2{pol) to W^ vpT (actually to W£ vpr in (F3p), we find for the 



longitudinally polarized beams, 



Pi(poi)P2(poi) W ^PT = 4hth 2 



p-q 



2 2 

y r 



sZ + 4 



p-q y r 



9l 



where we have introduced the variables, 

_ p-q 

y 



p-q 



(A.7) 



(A.J 



h-p h-q 

Hence the difference between the cross sections for the two-photon annhilation pro- 
cess with e + e~ beams polarized parallel and antiparallel to each other is given by 



Ax TT - da n 



d 1 d 1 2 



a" 



E[E' 2 <K\h-l 2 )pW 



4^2 2 2 

p-q y r 



+4 



Si 



f 

2/ r 



9l 



. (A.9) 



Especially for colliding beams, the laboratory is considered to be the cm. refer- 
ence frame. We have 

h = (E,0,0,E) , l 2 =(E,0,0,-E) , 

l[ = (E[, E^sinexcos^x, E^sine^in^x, E^cosOx) , (A.10) 
l 2 = (E' 2 , E' 2 sin92Cos(j)2, E' 2 sin92sm(j)2, —E 2 cos9 2 ) ■ 

where 9x, 4>x, and ir — 9 2 , (f) 2 are the polar and azimuthal angles for the final leptons 
l[ and 1' 2 , respectively. Then we obtain 



da n - da n 



a 



[(E + E[)(E + E' 2 ] 



E' 1 E 2 dE' 1 dE 2 dcos9 1 dcos9 2 d(f) 

7i E 2 p 2 q 2 (p ■ q) 

+(E + E[cos9x){E + E' 2 cos9 2 ) - E'xE' 2 sia9xsia9 2 cos(j) \g1 

+—E 2 E' l E2\{l - cos0i)(l - cos# 2 ) - 2sin0 1 sin0 2 cos0W 
p ■ q I J 



(AH) 

2 Our definition of W^ vpT and, therefore, gj and g\ , is such that they are proportional to 
e 2 (= 4ira), and not to e 4 in conformity to the nucleon case. 
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where = 0i — 02 , and 

q 2 = -2EE[(1 - cos0i) , p 2 = -2EE' 2 {\ - cosfl 2 ) , (A.12) 

p.q = [E - E[){E - E' 2 ) + (E -^cos0i)(£ - E' 2 cos6 2 ) - ^^sin^sin^cos^ . 



B Calculation of K^J 

In this appendix we present details of the calculation of K^J, the mixing anomalous 
dimension between the hadronic and photon operators, RT^ and -R" 3 ) 7 , for arbitrary 
n. The expressions of the operators RV 3 ^ and R%\ 7 are given in Eqs.( f4.H"D and 
( |3.13|) , respectively (we put aside the charge factor Q ch ). As a standard procedure, 
we introduce a light-like vector A M (A 2 = 0) to symmetrize the Lorentz indices and 
to eliminate the trace terms, and we define 

pn A — d <7 / 1 1' ,, /'ii-1 A A T? n A — D ff W'"/ 1 n-l A A fTZ ~l \ 

K {3)i ■ A = K (3)i A w • • • A Mn _ i; R {3h -/\ = K {3h A Ml • • • A /ln _ 1 . (B.1J 

Our first task is to evaluate the amputated two-point function with -R™ 3 ) 7 ■ A 
embedded between two photon fields at 0(1). We find 

< 7 |^3 )7 -A| 7 > = (0\T(A p (-p)R^ h -AA T (p))\0) Amp 

= i ^—^ {e^ppip ■ A)"- 1 - e pT ^A aPp p^p • A) n ~ 2 j (B.2) 

Next we calculate the one-loop diagram for the two-point function with R™ 3 y • A 
sandwiched by two photon fields, 

<7|i2fo, • A| 7 > = (0|T(A p (-p)i^ • AA T (p))\0) Amp , (B.3) 

which should be O(o). The operator R?fy is made up of four terms, VJ, Ui, VJ, and 
Ui, whose expressions are given in Eqs.( |4.15| - |4.18|) . Since VJ and U\ terms already 
have the QCD coupling constant g, their contributions are 0(g 2 a). So we work with 
Vi and Ui terms. The three-point "basic" vertices of V\ ■ A and U\ ■ A depicted in 
Fig. 8 are given, respectively, by 

V l>K = e l5 (k-A) l - 1 (k-A K -(k-A)g\)[(k-p)-A] n - 2 - 1 ^, (B.4) 
U hK = ie{k-A) l - l ^ K A a k p [{k-p)-A] n - 2 -%. (B.5) 
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The 0{a) contributions are coming from two diagrams shown in Fig.9. Inspecting 
the form of Vi )K we easily see that loop diagrams for V\ ■ A fail to produce a term 
which is proportional to (7|-R(3) 7 • A (7). In fact, after the loop integral we find that 
both Vi and terms give null result. 

On the other hand, one-loop diagrams for Ui and [7 n _i_; give terms proportional 
to (7|i?(3) 7 - AI7). The logarithmically divergent part of U\ contribution, for example, 
has the following form: 

i-< 7 |^ • A| 7 > = ^^ji-^m + 2,n- 0(7^ • A| 7 > InA 2 , (B.6) 

where Beta function B(l + 2, n — I) has appeared from the Feynman-parameter 
integral 

/ dxx l+1 (x- I)™- 1 - 1 = (-\) n - 1 - l B{l + 2,n-l) . (B.7) 
Jo 

Hence, adding together the U n ^\^i contribution, we find the mixing anomalous di- 
mension between R™ 3 y and Rfyj operators (apart from the quark charge factor) 

as 

KO) _ 24 



K 



n.l i 

n — 1 



l) n ~ l B(l + 2,n-l) + (-iy +1 B(n + 1-1,1 + 1) 



B.8 



In particular, for n = 3 and so I — 1, we have K^2 31 = —24-^ except for the quark 
charge factor. 

Now it is interesting to note the following identity: 

n-2 1 

y(n-l-l) x [(— l) n_ ' B(l + 2,n-l) + (-l) m B(n+l-l,l + l 

f=t n ~ 11 
1 l 

= (for odd integer n) , (B.9) 

n n+1 

which is a direct consequence of the relation (|4.21|) satisfied by the twist-3 operators. 
We know K^' m = from the null result of the mixing anomalous dimension between 
the mass operator R%\ m and -Rpw Also we know K^ E = since the photon matrix 
element of EOM operator R^ E vanishes. Thus we have a relation, 

n-2 

<l=E("-!-0<;- (B.10) 
1=1 
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The identity ( |B.9| ) assures that the relation indeed holds true. For n = 3 we have 



C Reanalysis of gl^fs 

In this appendix we reanalyze g1^n=3 , the third moment of ■ (The superscripts 

n = 3 and NS are omitted.) We choose R^ F ,R^i,R^ E as independent operators, 
replacing R(3) m with R(3)f- In these operator bases, the tree level coefficient func- 
tions are given by 

£ {3 )H1,0) = 1, £ (3 )i(1,0) = 0. (C.I) 

The relevant 2x2 anomalous dimension matrix 7^ has a form, 

/~(0) ~(o)\ 

z(p) - I "Iff 1f\ \ (r 9 \ 

7 - ~(o) ~(o) > 1^-4) 



v 



7!? 7!? 



with 



7® = 6C G -12C F , 7i? = 6C G (C.3) 

which are obtained from Eqs.(33)-(38) of Ref . ||47|| . (Note that we follow the conven- 
tion of Bardeen and Buras to define the anomalous dimension matrix which is the 
transposed one as given in Ref. ||47|| . ) Then, the eigenvalues of and corresponding 
projection operators, such that 7^ = \\Pi + A 2 P 2 , are found to be 

X 1 = 12C F , \ 2 = 6C G -^C F , (C.4) 



with 6 = Cf/[9(2Cf — Co))- Comparing the results in Eq.( |4.26| ), we note that 
Ai = 7^ vs ' and A 2 = ^i NS . This is the consequence of the fact that a different 
choice of the operator bases leads to different forms for the anomalous dimension 
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matrix but its eigenvalues remain the same. Now inserting these -E/(3)*(l, 0), X%, and 
Pi for i — F, 1 into the moment formula for ifl in Eq.( |4.9| ), we find that 



~1{NS) _ 
i/2,n=3 01 -2/ 



+ 







[- 




s 2 (6+l) 




U 2 








(£ 






U 2 



. 2 \ Ai/2/3 +l 



+ 



1 + 



2ft) 



4_ 



2 \ A 2 /2ft+l 



+ 



1 + ^2- 

1 ^ 2ft, 



2\ A 2 /2ft,+l 



(C.6) 



From Appendix B, we observe Kp^ = K± ^ for n = 3. Thus g^n=3 is reduced to 



_7(A r S) 1 #1 
<?2,n=3 ^ =2 



(0) 



y 1 + 2/3 



2 \ A 2 /2ft,+l 
2 



(C.7) 



which coincides with the expression for gln=z given in Eq.( (4.30| ) since A 2 = ^y^ 11 " . 



. (0)NS 



24 



References 

[1] A. V. Efremov and O. V. Teryaev, JINR Report NO. E2-88-287, Dubna, 1988; 
Phys. Lett. B240, 200 (1990). 

[2] S. D. Bass, Int. J. Mod. Phys. A7, 6039 (1992). 

[3] S. Narison, G. M. Shore and G. Veneziano, Nucl. Phys. B391, 69 (1993); 
G. M. Shore and G. Veneziano, Mod. Phys. Lett. A8, 373 (1993); 
G. M. Shore and G. Veneziano, Nucl. Phys. B381, 23 (1992). 

[4] A. Freund and L. M. Sehgal, Phys. Lett. B341, 90 (1994). 

[5] S. D. Bass, S. J. Brodsky and I. Schmidt, Phys. Lett. B437, 424 (1998). 

[6] D. P. Barber, in proceedings of the "Zeuthen Workshop on the Prospects of Spin 
Physics at HERA", DESY 95-200, p.76, eds. J. Bliimlein and W. D. Nowak. 

[7] M. Stratmann and W. Vogelsang, Z. Phys. C74, 641 (1997). 

[8] R. W. Brown and I. J. Muzinich, Phys. Rev. D4, 1496 (1971). 

[9] M. A. Ahmed and G. G. Ross, Phys. Lett. B59, 369 (1975). 

[10] K. Sasaki, Phys. Rev. D22, 2143 (1980); Prog. Theor. Phys. Suppl. 77, 197 
(1983). 

[11] A. V. Manohar, Phys. Lett. B219, 357 (1989). 

[12] M. Stratmann and W. Vogelsang, Phys. Lett. B386, 370 (1996). 

[13] M. Gliick, E. Reya and C. Sieg, Phys. Lett. B503, 285 (2001); Eur. Phys. J. 
C20, 271 (2001). 

[14] W. A. Bardeen and A. J. Buras, Phys. Rev. D20, 166 (1979). 

[15] D. W. Duke and J. F. Owens, Phys. Rev D22, 2280 (1980). 

[16] T. Uematsu and T. F. Walsh, Phys. Lett.lOlB, 263 (1981). 

25 



[17] T. Uematsu and T. F. Walsh, Nucl. Phys. B199, 93 (1982). 
[18] G. Rossi, Phys. Rev. D29, 852 (1984). 

[19] M. Drees and R. M. Godbole, Phys. Rev. D50, 3124 (1994). 

[20] M. Gliick, E. Reya and M. Stratmann, Phys. Rev. D51, 3220 (1995); 
Phys. Rev. D54, 5515 (1996); 

M. Gliick, E. Reya and I. Schienbein, Phys. Rev. D60, 54019 (1999). 

[21] B. L. Ioffe and A. Oganesian, Z. Phys. C69, 119 (1995). 

[22] G. A. Schuler and T. Sjostrand, Phys. Lett. B376, 193 (1996). 

[23] P. Mathews and V. Ravindran, Int. J. Mod. Phys All, 2783 (1996) . 

[24] K. Sasaki and T. Uematsu, Phys. Rev. D59,114011 (1999) ; Phys. Lett. B473, 
309 (2000); Eur. Phys. J. C20, 283 (2001). 

[25] E143 Collaboration, K. Abe et al, Phys. Rev. Lett. 76, 587 (1996). 

[26] E155 Collaboration, P. L. Anthony et al, Phys. Lett. B458, 529 (1999). 

[27] S. Wandzura and F. Wilczek, Phys. Lett. B 172, 195 (1977). 

[28] A. Gabrieli and G. Ridolfi, Phys. Lett. B 417, 369 (1998). 

[29] H. Burkhardt and W. N. Cottingham, Ann. Phys. 56, 453 (1970). 

[30] W. A. Bardeen, A. J. Buras, D. W. Duke, and T. Muta, Phys. Rev. D18, 3998 

(1978) . 

[31] J. Kodaira, S. Matsuda, K. Sasaki and T. Uematsu, Nucl. Phys. B159, 99 

(1979) . 

[32] E. V. Shuryak and A. I. Vainshtein, Nucl. Phys. B199, 451 (1982); B201, 141 
(1982). 

[33] A.P Bukhvostov, E.A. Kuraev and L.N. Lipatov, Sov. Phy. JETP 60, 22 (1984). 

26 



[34] P.G. Ratcliffe, Nucl.Phys. B264, 493 (1986). 

[35] I.I. Balitsky and V.M. Braun, Nucl. Phys. B311, 541 (1988/89). 

[36] X. Ji and C. Chou, Phys. Rev. D42, 3637 (1990). 

[37] R. L. Jaffe, Comments Nucl. Part. Phys. 19, 239 (1990). 

[38] R. L. Jaffe and X. Ji, Phys.Rev. D43, 724 (1991). 

[39] J. Kodaira, T. Uematsu and Y. Yasui, Phys. Lett. B344, 348 (1995). 

[40] J. Kodaira, K. Tanaka, T. Uematsu and Y. Yasui, Phys. Lett. B387, 855 (1996). 

[41] D. Miiller, Phys. Lett. B407, 314 (1997). 

[42] V. M. Braun, G. P. Korchemsky and A. N. Manashov, Phys.Lett. B476, 455 
(2000); Nucl. Phys. B597, 370 (2001); B603, 69 (2001). 

[43] A. V. Belitsky, X. Ji, W. Lu, and J. Osborne, Phys. Rev. D63, 094012 (2001). 

[44] J. Kodaira, T. Nasuno, H. Tochimura, K. Tanaka, and Y. Yasui, Prog. Theor. 
Phys. 99, 315 (1998); |hep-ph/ 97 10451 , in Proceedings of 2nd Topical Workshop 



on Deep Inelastic Scatterings off Polarized Targets (1997) Zeuthen, p. 210. 

[45] A. Ali, V. M. Braun and G. Hiller, Phys. Lett. B266, 117 (1991). 

[46] V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 675 (1972); L. N. Li- 
patov, ibid. 20, 94 (1975); Y. L. Dokshitser, Sov. Phys. JETP 46, 641 (1977); 
G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977). 

[47] K. Sasaki, Phys. Rev. D58, 094007 (1998). 

[48] V. M. Budnev, I. F. Ginzburg, G. V. Meledin and V. G. Serbo, Phys. Rep. 15 
(1975) 181. 



27 



Figure Captions 



Fig. 1 Deep inelastic scattering on a polarized virtual photon in polarized e + e~ 
collision, e + e~ — > e + e~+ hadrons (quarks and gluons). The arrows indicate 
the polarizations of the e + and e~. The mass squared of the "probe" ("target") 
photon is -Q 2 (-P 2 ) (A 2 < P 2 < Q 2 ). 

Fig. 2 Forward scattering of a virtual photon with momentum q and another virtual 
photon with momentum p. The Lorentz indices are denoted by /i, v, p, r. 

Fig. 3 The Box-diagrams contributing to gj and g\ in the pure QED interaction. 

Fig. 4 The Box-diagram contributions to g1{x, Q 2 , P 2 ) (dashed line), g1{x, Q 2 , P 2 ) 
(solid line) and gl(x, Q 2 , P 2 ) (dash-2dotted line) for Q 2 = 30 GeV 2 and P 2 = 1 
GeV 2 for Nj = 3. The 2x — 1 line shows the leading logarithmic term of gj. 

Fig. 5 Triangle diagrams which contribute to the anomalous dimension describing 
the mixing between the twist-3 quark-bilinear operator R?$\ q and the photonic 
operator R? 3 \ 7 - 

Fig. 6 The third (n = 3) moment of g]>(x, Q 2 , P 2 ) in units of a/ff as a function of 
q2 f or p2 _ j Q e y2 wit j 1 jy^ _ 3^ rj^g dash-dotted (dashed) line corresponds 

to the flavor nonsinglet (singlet) component. The solid line represents the sum 
of the two components, the total gl n=3 /(a/ir). 

Fig. 7 The Box-diagram (dashed line) and the QCD LO (solid line) contributions 
for large N c to the flavor nonsinglet photon structure function g1^ NS ^ (x, Q 2 , P 2 ) 
for Q 2 = 30 GeV 2 and P 2 = 1 GeV 2 for N f = 3. 

Fig. 8 The tree-level three-point vertices of V\ and U\. 

Fig. 9 One-loop diagrams at O[o) contributing to the Green's functions of V\ and 
Ui with two photons as external lines. 
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